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Abstract. We introduce a notion of chain of evolution algebras. The 
sequence of matrices of the structural constants for this chain of evo- 
lution algebras satisfies an analogue of Chapman-Kolmogorov equation. 
We give several examples (time homogenous, time non-homogenous, pe- 
riodic, etc.) of such chains. For a periodic chain of evolution algebras 
we construct a continuum set of non-isomorphic evolution algebras and 
show that the corresponding discrete time chain of evolution algebras 
is dense in the set. We obtain a criteria for an evolution algebra to be 
baric and give a concept of a property transition. For several chains of 
evolution algebras we describe the behavior of the baric property de- 
pending on the time. For a chain of evolution algebras given by the 
matrix of a two-state evolution we define a baric property controller 
function and under some conditions on this controller we prove that the 
chain is not baric almost surely (with respect to Lebesgue measure). We 
also construct examples of the almost surely baric chains of evolution 
algebras. We show that there are chains of evolution algebras such that 
if it has a unique (resp. infinitely many) absolute nilpotent element at a 
fixed time, then it has unique (resp. infinitely many) absolute nilpotent 
element any time; also there are chains of evolution algebras which have 
not such property. For an example of two dimensional chain of evolution 
algebras we give the full set of idempotent elements and show that for 
some values of parameters the number of idempotent elements does not 
depend on time, but for other values of parameters there is a critical 
time t c such that the chain has only two idempotent elements if time 
t > t c and it has four idempotent elements if time t < t c . 
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1. Introduction 

In this paper we consider some classes of non-associative algebras. There 
exist several classes of non-associative algebras (baric, evolution, Bernstein, 
train, stochastic, etc.), whose investigation has provided a number of signifi- 
cant contributions to theoretical population genetics. Such classes have been 
defined at different times by several authors, and all the algebras belonging 
to these classes are generally called "genetic". Etherington introduced the 
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formal language of abstract algebra to study of genetics in his series of sem- 
inal papers [2]-|4j . In recent years many authors have tried to investigate the 
difficult problem of classification of these algebras. The most comprehensive 
references for the mathematical research done in this area are [H ]l l2 | IT5 l ll6| . 

In [11] an evolution algebra A associated to the free population is intro- 
duced and using this non-associative algebra many results are obtained in 
explicit form, e.g. the explicit description of stationary quadratic operators, 
and the explicit solutions of a nonlinear evolutionary equation in the absence 
of selection, as well as general theorems on convergence to equilibrium in 
the presence of selection. 

In |15j a new type of evolution algebra is introduced. This evolution 
algebra is defined as follows. Let (E, •) be an algebra over a field K. If it 
admits a basis e±, e^i ■ ■ ■ > such that e, • ej = 0, if i ^ j and • ej = a^e^, 
for any i, then this algebra is called an evolution algebra. 

In this paper by the term evolution algebra we will understand a finite 
dimensional evolution algebra E (as mentioned above) over field R. 

Evolution algebras have the following elementary properties (see [15]): 
Evolution algebras are not associative, in general; they are commutative, 
flexible, but not power-associative, in general; direct sums of evolution alge- 
bras are also evolution algebras; Kronecker products of evolutions algebras 
are also evolution algebras. 

The concept of evolution algebras lies between algebras and dynamical 
systems. Algebraically, evolution algebras are non- associative Banach al- 
gebra; dynamically, they represent discrete dynamical systems. Evolution 
algebras have many connections with other mathematical fields including 
graph theory, group theory, stochastic processes, mathematical physics, etc. 

In the book [15], the foundation of evolution algebra theory and appli- 
cations in non-Mendelian genetics and Markov chains are developed, with 
pointers to some further research topics. 

In [13] the algebraic structures of function spaces defined by graphs and 
state spaces equipped with Gibbs measures by associating evolution algebras 
are studied. Results of [13] also allow a natural introduction of thermody- 
namics in studying of several systems of biology, physics and mathematics 
by theory of evolution algebras. 

The paper is organized as follows. In Section we give main definitions 
related to a chain of evolution algebras. Therein we give several examples 
(time homogenous, time non-homogenous, periodic, etc.) of such chains. 
For a periodic chain of evolution algebras we construct a continuum set 
of non-isomorphic evolution algebras and show that the corresponding dis- 
crete time chain of evolution algebras is dense in the set. In Section 3 we 
obtain a criteria for an evolution algebra to be baric. The concept of a 
property transition is introduced in Section 4. This section also contains 
several chains of evolution algebras for which we describe the behavior of 
the baric property depending on the time. For a chain of evolution algebras 
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given by the matrix of a two-state evolution we define a baric property con- 
troller function and under some conditions on this controller we prove that 
the chain is not baric almost surely (with respect to Lebesgue measure). 
We also construct examples of the almost surely baric chains of evolution 
algebras. We show that there are chains of evolution algebras such that if 
it has a unique (resp. infinitely many) absolute nilpotent element at a fixed 
time, then it has unique (resp. infinitely many) absolute nilpotent element 
any time; also there are chains of evolution algebras which have not such 
property. In the last subsection for an example of two dimensional chain of 
evolution algebras we give the full set of idempotent elements and show that 
for some values of parameters the number of idempotent elements does not 
depend on time, but for other values of parameters there is a critical time 
t c such that the chain has only two idempotent elements if time t > t c and 
it has four idempotent elements if time t < t c . 

2. Definition and examples of CEA 

Consider a family jl?[ s >*] : s, t 6 K, < s < i} of n-dimensional evolution 
algebras over the field M, with basis ci, . . . , e n and multiplication table 

n 

(2.1) aei = M\ s,£[ = ^ a tj i = 1, ■ ■ ■ , n; aej = 0, i ^ j. 

i=l 

Here parameters s,t are considered as time. 

Denote by M\ s ^ = (at*.'*') -the matrix of structural constants. 

V / i,j=l,...,n 

Definition 2.1. A family ji£[ s >*] : s, t G R, < s < i} of n-dimensional 
evolution algebras over the field M. is called a chain of evolution algebras 
(CEA) if the matrix M\ s,t ^ of structural constants satisfies the Chapman- 
Kolmogorov equation 

(2.2) M [s ' t] = M [s ' t] M [tA , for any s<T<t. 

If pi is a projection map of E^ s,t \ which maps every element of E^ s,t ^ to 
its ei component, then equation (12. 2p can be written as 

n 

(2.3) Mf A =J2Pj( M l S ' T] Wj r ' t \ for an Y s<r<t. 

i=i 

Definition 2.2. A CEA is called a time-homogenous CEA if the matrix 
Ji4^ s '^ depends only ont — s. In this case we write A4\ l ~ s \ 

Definition 2.3. A CEA is called periodic if its matrix M\ s,t * is periodic with 
respect to at least one of the variables s, t, i.e. (periodicity with respect to 
t) M^ s ' t+P ^ = M\ sA for all values of t. The constant P is called the period, 
and is required to be nonzero. 
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Remark 2.4. In general, an algebra A^'^ can be given by a cubic matrix 
M^'^ = (a&u) of structural constants. Our Definition can be 

extended to A^'^ using analogues of the Chapman- Kolmogorov equations for 
quadratic operators (see ®Yfy\lJ$). Since in the general case there are two 
types of the Chapman- Kolmogorov equations: type A and type B JB^, one 
also can define two types of chain of (general) algebras using the Chapman- 
Kolmogorov equations of type A and type B, respectively. In this paper we 
shall only consider CEA, which is more simple than general case, because it 
is defined by quadratic matrices. 

The CEA corresponding to a Markov process. 

Let {.A/f ' s '*l, < s < t} be a family of stochastic matrices which satisfies 
the equation (|2.2|) . then it defines a Markov process. Thus we have 

Theorem 2.5. For each Markov process, there is a CEA whose structural 
constants are transition probabilities of the process, and whose generator set 
(basis ) is the state space of the Markov process. 

If Ai^'^ does not depend on time (i.e. = A4) then the CEA contains only 
one evolution algebra E. Note that for a Markov chain defined by Ai the 
corresponding E has been studied in [15] . 

Now we shall give several concrete examples of CEA. 

Example 1. To show a time dependent CEA we use the following example 
of time homogenous Markov process (see [10J) : for n = 3 consider 

4? = ^e-t A *cos(ai) + i z = 1,2,3; 
41 = 4! = 41 = e ~^* ( -Jg sin(at) - ^ cos(at)) + ^; 

41 = 41 = 41 = -e~'^ At sin(art) + i cos(ai)^ + i, 

where A > 0, a = ^-A. 

Let E^\t > be the corresponding CEA. It is easy to see that has 
an oscillation behavior depending on time t. Moreover limt_j. +00 = E, 
where E is an evolution algebra with the multiplication table 

e\ = el = ef = -(ei + e 2 + e 3 ), e i e j = 0,i ± j. 

The CEA corresponding to a family of matrices which do not define a 
process. 

Example 2. We shall give a time homogenous CEA which are different 
from CEAs corresponding to Markov processes. For n = 2 take 

„[*] __[*] __[*]. _M __M _ b [*] 

"11 — "22 ~~ u J u 12 — a 21 ~~ u ■ 

Then equation (|2.2p is equivalent to 

a [t)= a \T) a [t-T] +b [T} h \t-T). 
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b [t] =a [r] b [t-r] +b lr] a [t-r]_ 

Denote f(t) = + &M, <p(t) = — then the last system of functional 
equations can be written as 

f(t) = f(r)f(t-r), <p{t) = <p{T)<p{t-T). 

Both these equations are known as exponential Cauchy equation and the 
system of equations has solution fit) = A', ip(t) = fi l , where A,/x > 0. 
Consequently, = ^(A* + /J), feM = ^(A* — /U*). But this solution does not 
define any Markov process, in general. 

Let E^\t > be the corresponding CEA. Depending on parameters A 
and fi we get distinct behavior of E^ for t — > +oo, i.e. we have 

( E if < X,n < 1, 
Ex if A = n = 1, 
E 1/2 if A = 1,0 < n < 1, 
£_ 1/2 if n = 1,0 < A < 1, 
E'oo otherwise, 



lim #M 



where E'o is an evolution algebra with zero multiplication; E\ is an evolution 
algebra with multiplication table 

el = e 1 ,el = e 2 , eie 2 = 0; 
E 1 / 2 is an evolution algebra with multiplication table 

e i = e 2 = ^( e i + e 2), e 1 e 2 = 0; 
E_i/ 2 is an evolution algebra with multiplication table 

e? = ^(ei-e 2 ), el = --(ei - e 2 ), eie 2 = 0; 

and E'oo is a vector space which has "infinity multiplication" , or we can say 
that in E^ an algebra structure is not defined. This example shows that a 
limit of a CEA can be non evolution algebra. 

Example 3. A two-state evolution. Now we shall give an example of 
time non-homogeneous CEA, the matrix of structural constants of which 
also does not define any (time non homogenous) Markov process in general. 

Consider n = 2 and matrix M\ s '^ = (af/^) with 

V V Ji,j=l,2 

a[f = \ (1 + a{s, t) + P{s, t)) , a[f = \ (1 - a(«, t) - (3(s, *)) , 
(2 ' 4) 4i* ] = 5 (! + a M) " , a!** 1 = i(l-a(M) +£(*,*))■ 
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In this case the equation (12. 2|) is equivalent to (see [9]) 

a(s, t) = a(r, t) + a(s, r)/3(r, t), 
( 2-5 ) P(s, t) = P(s, r)/3(r, t), s<T<t. 

The second equation of the system (|2,5p is known as Cantor's second equa- 
tion, it has very rich family of solutions: f3(s,t) = |^y, where $ is an 
arbitrary function with <&(s) 7^ 0. Using this function (3 for the function a 
we obtain 

a(s,t) a(r, i) a(s,r) 
= $(t) + $(r) ' 

Now denote 7(s,i) = ifey then the last equation gets the following form 

j{s,i) = 7(s,r) +7(r,t). 

This equation is known as Cantor's first equation which also has very rich 
family of solutions: j(s,t) = ^(t) — ^(s), where ^ is an arbitrary function. 

Hence a solution M\ s $ = (aj*-'*' ) to the equation (|2,2p is given by 



i,J=l,2 



•^-ifi-<6(t)(*(t) -*(,)) 



'12 



2 V 

».*]- 1 +< & (t)(*(t)_ ¥( a )) Ml 



'21 



2 V *(«) 



Let ^[ s '*], < s < t be the corresponding to this solution CEA. This CEA 
varies by two parameters, for example, if t = s we get E^ = E with 
multiplication table ef = ei, e| = e%, e\e 2 = 0. Moreover, choosing functions 

and one can variate the limit behavior of the CEA. For example, if 
and such that lim^+oo <3>(t)^(t) = limi_j. +00 <I>(i) = 0, then for a fixed 
s we have lim^+oo E^- S,t ^ = Ei/ 2 , where Ei/ 2 is an evolution algebra with 
multiplication table 

e \ = e 2 = \i e i + e 2)i eie 2 = 0. 

Example 4. A n- dimensional time non-homogenous CEA. Here for ar- 
bitrary n we shall give an example of time non-homogenous CEA. Let 
{A^\ t > 0} be a family of invertible (for all t), n x n matrices. Define 
the following matrix 

M [s,t) = aW(A®)~ 1 , 
where (A^) -1 is the inverse of A^. 



A CHAIN OF EVOLUTION ALGEBRAS 



7 



This matrix satisfies the equation (|2.2p . Indeed, using associativity of the 
multiplication of matrices we get 



] M [r,t] = A [s] ^(^M)-i^M^ (yiM)- 1 = A^iA®)- 1 = M [s ' t] . 



Thus each family (with one parameter) of invertible n x n matrices defines 
a CEA I?[ s '*] which is time non-homogenous, in general. But will be a time 
homogenous CEA, for example, if A$ is equal to tth power of an invertible 
matrix A. 

Construction of a family of invertible n x n matrices A® is not difficult, 
for example, one can take A® triangular n x n matrix of the form 



A lt] 



11 







a 



M „[*] 



21 



'22 








'31 



'32 



,[*] 



\ 








a [t] 

"■7171. 



b n2 ■ ■ ■ "nn-1 ""Tin ) 

which is called lower triangular matrix or one can take an upper triangular 
matrix. Then the matrices are invertible iff a$ ^ 0, for all i = 1, . . . , n and 
t. So this example also gives a very rich class of CEAs. 

Example 5. Periodic CEA. To get a periodic CEA, we can consider 
the constructed in Example 3, and choose $ and \& as periodic (non- 
constant) functions. Then corresponding CEA is periodic. In this case for 
any fixed s, the limit lim^+oo E^'^ does not exist in general, moreover its 
set of limit points (evolution algebras) can be a continuum set. We shall 
make this point clear as follows. Construct a time homogenous CEA which 
is periodic. Consider n = 2 take 



i [t] 
hi 



i [t] 



,!*]■ 



'12 



'21 



Then equation (|2.2p is equivalent to 

a [*] =a M [*-r]_ 6 [r] 6 tt-r]. 
6 [*] =0 W 6 [t-r] +6 [r) a [t-r]_ 

This system reminds the following identities 

cos t = cos r cos(t — r) — sin r sin(t — r) 
sin t = cos r sin(t — r) + sin r cos(t — r) 

Consequently, one solution = ( a^- ) t o equal ion (2.2) is 

(2.6) M [t] = 



[t}\ 

\ ij Ji,j=l,2 

cos t sin t 
- sin t cos i 
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Since that matrix is periodic with period P = 2ir, the corresponding CEA 
E$ is also periodic. Moreover this CEA is very interesting: for arbitrary 
2-dimensional evolution algebra E+ , or E~ , a G [—1,1] with structural 
constants matrix 



Ml 



± ( a ±\/l — a, 



=FVl - a 



2 



a 



respectively, there is a sequence t n = t n (a) of times such that lim^oo E^ = 
E+ or E~ . We have ^ E^ if a 7^ b. Moreover the following is true 

Proposition 2.6. 1) For any a, b G [—1,1], a 7^ ±6, i/ie algebras E+ and 
E^ are not isomorphic. The algebras E^ and E^ a are isomorphic. 

2) For any a,b £ [—1, 1], a / ±b, the algebras E~ and E~ are not iso- 
morphic. The algebras E~ and EZ a are isomorphic. 

/ a p \ 

Proof. 1) Let ip = I ^ I be an isomorphism of the evolution algebra 

E+ to the evolution algebra E^ . Here det(y) 7= 0. By the multiplication ta- 
ble of the evolution algebras, we get the following relation between matrices 
M+ and M+: 

b det(p) 

(aa + v 7 ! — a 2 /3)/3 2 — (a/3 — v 7 ! 



(a<5 - v / T 3 ^7)7 2 - (07 + v 7 !^ 2 "^ 2 (aa + v 7 ! 3 ^ 2 ^ 2 - (a/3 - vT 

Since det(7W+) = 1, it is easy to see that there are two classes of isomor- 
phisms: 



For the class Ci the matrix A4^ must satisfy the following 

/ /, w'T - la \ / v 1 - 

Ml 



-VT^b 2 b 



-VT^ 2 ^ a5 



From this equality we get a = 5 = y = ^ if a 7^ 0, ±1 which is 
satisfied iff a = 6. Hence the isomorphisms from the class C\ can not give 
an isomorphism from A4+ to M£ ■ For a = we get 6 = 0. One can take 
a = 5 = =Fl if a = ±1 and b = =f1. Hence E^ is isomorph to E^- v 
For the class C2 the matrix M^ must satisfy the following 

. . . / „fl -/I 12/3 2 



A CHAIN OF EVOLUTION ALGEBRAS 







From this equality we get f3 = 7 = — y ]__^ = f if a 7^ 0, ±1 which is 
satisfied iff a = —b. Hence the isomorphisms from the class C2 can only give 
an isomorphism from Ai^ to Al_ a . 

2) The proof of 2) is similar to the proof of 1). □ 

Consider now discrete time n, n E N and the CEA {.E^l, n E N} given by 
matrix (|2.6p . 

Proposition 2.7. T/ie discrete time CEA E^ n \ n E N, is dense in the 
set {E^, a E [—1,1]} of evolution algebras, i.e. for an arbitrary evolution 
algebra E^ there exists a sequence {nk}k=i,2,... of natural numbers such that 
]im k ^ 00 EM=E+ orE~. 

Proof. It is known that the sequences {sin n} and {cos n}, n E N, are dense in 
[— 1, 1] (see e.g. [5]). Hence for any a E [—1, 1] there is a sequence {nk\k=i,2,... 
of natural numbers such that lim^oo cos(n^) = a. The same sequence can 
be used to get lim^oo E^ = E+ or E~ . □ 



3. A CRITERION FOR AN EVOLUTION ALGEBRA TO BE BARIC 

A character for an algebra A is a nonzero multiplicative linear form on A, 
that is, a nonzero algebra homomorphism from A to R |llj . Not every alge- 
bra admits a character. For example, an algebra with the zero multiplication 
has no character. 

Definition 3.1. A pair (^4,0") consisting of an algebra A and a character 
a on A is called a baric algebra. The homomorphism a is called the weight 
(or baric) function of A and o~(x) the weight (baric value) of x. 

In [11] for the evolution algebra of a free population it is proven that 
there is a character o~(x) = ^i 3 ^' therefore that algebra is baric. But the 
evolution algebra E introduced in [TB] is not baric, in general. The following 
theorem gives a criterion for an evolution algebra E to be baric. 

Theorem 3.2. An n-dimensional evolution algebra E, over field the M, 
is baric if and only if there is a column (au , . . . , a n i ) T of its structural 
constants matrix M. = (ay)i 3 =i n , such that Oj j 7^ and au = 0, for all 
i 7^ iQ. Moreover, the corresponding weight function is o~(x) = ai Q i Q Xi . 

Proof. Necessity. Take x,y E E with x = Y^i=l x * e *' V = Y17=i V^i- Assume 
o{x) = Yli=i a i x i-> x £ E is & character. We have 

n I n \ ran 

<r{vv) = ^2 ^2 a ij a j x iVu o-(x)a(y) = ^ ^ a^ x^-. 
i=i \j=i J i=i j=i 

From a(xy) = o(x)a{y) we get 

(3.1) aiCtj = for any i ^ j, i, j = 1, . . . , n; 
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n 

(3.2) a ij a j = a l f° r an Y i = I, ■ ■ ■ ,n. 

3=1 

It is easy to see that the system (|3. 1 j) has a solution a = («i, . . . , a n ) with 
af + • • • + > if and only if exactly one coordinate of a, say aj , is not 
zero, and all others are zeros. Substituting this solution in (|3,2p we get 

au ai = 0, if i ± io, i = 1, . . . , n; 

ai io«io = «i > if * = *o- 
From the last equations we get aj j / 0, ajj = 0, for all i ^ io and 
Oti = a io«o- 

Sufficiency. Assume there is a column (au , . . . , a n i ) T ', such that aj j 7^ 
and ajj = 0, for all z 7^ io- Then it is easy to see that a(x) = ai i Xi is a 
weight function, therefore E 1 is a baric evolution algebra. □ 

A baric algebra A may have several weight functions. As a corollary of 
Theorem 13.21 we have 

Corollary 3.3. If the matrix Ai, mentioned in Theorem \3.2l has several 
columns (au j , . . . , a^) , j = i\, . . . , i m , m < n, which satisfy conditions 
of Theorem \3.2\ then the evolution algebra E has exactly m weight functions 

0~(X) — fHjij^ij} 3 — *lj • • • lira' 

There are two types of trivial evolution algebras [15J: zero evolution alge- 
bra, which satisfies = for all i, j = 1, . . . ,n; non-zero trivial evolution 
algebra, which satisfies ejej = for all i ^ j and ef = auei, where an G E is 
non-zero for some i = 1, . . . , n. By Theorem 13.21 we conclude that the zero 
evolution algebra is not baric, but any non-zero trivial evolution algebra is 
a baric algebra. Moreover, there are baric evolution algebras which are not 
trivial. 

4. Property transition 

If a system has parameters (as usually like: temperature, time, interac- 
tion, etc.) then a property of the system can variate by a parameter. For 
example, the behavior of phases (states) of a system in physics, depends on 
temperature T > 0, if for some values of T there is a unique phase and for 
other values there are several phases, then the physical system has a phase 
transition [8]. Similar transitions of a property can be seen for systems of 
biology, chemistry, etc. Here we shall define a notion of property transition 
for CEA. 

Definition 4.1. Assume a CEA, E^ s,t \ has a property, say P, at pair of 
times (so,to); we say that the CEA has P property transition if there is a 
pair (s,t) 7^ (sQ,to) at which the CEA has no the property P. 



Denote 



T={(s,t) :0<s<t}; 
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T P = {{s,t) E T : E [s ' t] has property P}; 

Tp 1 = T\T P = {(s,t) £T : E [s ' t] has no property P}. 

Definition 4.2. We call the set 
Tp-the duration of the property P; 
Tp-the lost duration of the property P; 

The partition {Tp,Tp} of the set T is called P property diagram. 

For example, if P =commutativity then since any evolution algebra is 
commutative, we conclude that any CEA has not commutativity property 
transition. 

4.3. Baric property transition. Since a CEA is not a baric algebra, in 
general, using Theorem 13.21 we can give baric property diagram. Let us do 
this for the above given Examples 1-4. 

Example 1'. For the case of Example 1, by Theorem 13.21 we have that 
is baric iff 

oSJ = |e-i A ooB(at) + i = l. 

This has unique solution t = 0. Consequently, Tbaric = {0}, 7^ aiic = {t : 
t > 0}. Thus the CEA E^ is baric (even non zero trivial) evolution algebra 
only at initial time, and it loses baricity as soon as the time turned on. 
Example 2'. In Example 2, using Theorem 13.21 we obtain that 



baric 



{0} if A^/i; 
T if A = /i. 

Thus the CEA has not baric property transition if A = fi, and it has a 
baric property transition, as in Example 1', if A 7^ [i. 

Example 3'. Baric property transition for a two-state evolution. Since in 
case of Example 3, we have a rich class of CEA here we shall give a special 
theory of the baric property transition. Using Theorem 13.21 we obtain that 
Tbaric is the set of (s,t) such that 

1 + mm) - *oo) - 1|| = or 1 - mm) - *w) - 1|| = o. 

These equations can be rewritten as 

9(t) = e( s ), e-(t) = e-( s ), 

where 

(4.1) ^ ) = _L_ + m ,- (t ) = _L__^). 

Thus 

"Tbaric "Tbaric \™) U Tbaric (fi )> 

here T baric (0) = {(s,t) £ T : 0(t) = 9(s)}. 
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Remark 4.4. To describe the set 7baric one has to describe the sets 7baric(#) 
and 7b a ric(0 - ); both of which are defined by the parameter functions $ and 
. Note that if we replace $ with — or \& with then these sets transfer 
to each other. Since $ and are arbitrary functions, it will be enough to 
describe only 7baric(#) f or arbitrary 9. Thus in the sequel of this subsection 
we shall deal with description of Tbaric (#)■ 

The function 9{t) is called baric property controller of the CEA. Because, 
it really controls the baric duration set, for example, if 9 is a strong monotone 
function then the duration is "minimal", i.e. the line s = t, but if 9 is a 
constant function then the baric duration set is "maximal" , i.e. it is T ■ Since 
$ and are arbitrary functions, we have a rich class of controller functions, 
therefore we have a "powerful" control on the property to be baric. 

For a special choose of 9 we have 

Proposition 4.5. // $(t) = \ l , A > and ^(t) = ct, c£ R. Then 

7baric(#) = 7Laric(A,c) = {(s,t) : S = t}U 

' if < A < 1, c> In A; or 

A > 1, c G (-oo,0] U [In A, +oo), 
* {{s, t) : < s < t c , t c < t < t' c , 6(s) = 9(t)} if < A < 1, c < In A; or 

A > 1, c G (0, In A), 

where t c and t' c serve as critical times, which defined by t c = hi ) 
and t' c > is a unique solution to 9(t' c ) = 1. 

Proof. Under the conditions of the proposition we have 9{t) = A _< + ct, and 
the simple analysis of the equation 9(s) = 9{t) for this 6 gives the full set 

7baric(A,c). □ 

In Figure [U the baric property diagram is given 




Figure 1. The baric property diagram for 9{t) = A t + ct. 
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As a corollary of Proposition 14.51 we have 

Corollary 4.6. 1) For any fixed s, with < s < t c (resp. t c < s < t), the 

time t has two (resp. one) critical values: t^ = s (resp. s) and t^ which 

(2) 

is a unique solution of9(t y c >) = 0{s). 

2) For any fixed t, with < s < t < t c or t' c < t (resp. t c < t < t' c ), the 

time s has one (resp. two) critical values: Sc = t (resp. Sc = t and s^f 1 

(2) 

which is a unique solution of6(t K c >) = 6(s)) . 

Let us discuss some more examples of the controller 9. If 9{t) = tan(i) 
then tan(s) = tan(i) has solution t = s + trk, k £ Z. The intersection of this 
family of lines with T gives the family of half lines, i.e. 



If 9{t) = sin(t) then sin(s) = sin(t) has two family of solutions: s = t + 2nk, 
k G Z and s = — t + (2k + l)7r, k € Z. The intersection of these families of 
lines with T is 



T b aric(sm(t)) = (J {(s,t)eT :t = s + 2irk or t = -s + (2k + l)vr} . 



In all above considered examples we obtained a set 7b a ric(6) which has 
zero Lebesgue measure. But there is controllers for which this set has non- 
zero Lebesgue measure, for example, if 9{t) is a controller function with the 
graph as shown in Figure [21 then the corresponding baric property diagram 
is as shown in Figure El Thus any "constant part" of the graph of the 
controller gives a full triangle in the diagram, moreover, any "non-constant 
part" gives several curves. In this case the set % a ric(6) has a non-zero 
Lebesgue measure. 

Let 9(t) = D{t) be the Dirichlet function defined by 



In this case we have a rich set of baric property duration, i.e. 

%aric{D(t)) = {(s,t) G T- t and s rational} U {(s, t) E T : t and s irrational}. 

Definition 4.7. A function 9 defined on R is called a function of countable 
variation if it has the following properties: 

1. it is continuous except at most on a countable set, (which is denoted 
by X c = {x\,X2, ■ ■ ■}), it has only jump-type discontinuities (denote the 
one-sided limit from the negative direction by 9(x^) and from the positive 
direction by 9(xf), i = 1, 2, ... ); 

2. it has at most a countable set of singular (extremum) points (which is 
denoted by X e = {y 1 ,y 2 , . . . }). 



7^ ar j c (tan(t)) = (J {(s, t) £ T : s = t - irk} . 



fc=0,l,2,... 



fc=0,l,2,... 




1 if t rational; 
if t irrational. 
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S 




t\ t 2 t$ titfrts tytgtgtiQ tu t 



Figure 2. An example of controller 6. 




ti t- 2 t$ ti t 5 tfi ty t$ tg t w tu t 



Figure 3. The baric property diagram for the controller 9 
with graph as in Figure [5J 



Note that any function of countable variation has not "constant parts 
its graph. 
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The following theorem gives a characteristics of the baric property dura- 
tion set. 

Theorem 4.8. If the controller 9 (see is a function of countable 

variation, then the baric duration set 7baric(#) has zero Lebesgue measure, 
that is the corresponding CEA is not baric almost surely. 

Proof. Using the (finite or infinite) sequences X c and X e we construct the se- 
quences {t^.}fc=i,2,..., i = 1, 2, . . . with 0(tr k ) = 0(x~) for all k; {i+ 9 } ? =i,2,..., 
i = 1,2,.'.. with 0(i+) = e{x+) and {^} i= i, 2 ,..., j = 1,2,.'.., where 
d(tjj) = 0(yj) for all I. Now define a sequence {ti}j=i,2,..., with t\ < t<i < 
t% < . . . as follows 

{ti}i=l A ... =X c UX e \J ({/,,}, ,.•>... U {i+ },=i, 2 ,...) U{^,/}i=lA.... 

Since 8 is a function of countable variation, the sequence {ii}i=i,2,... is at 
most countable. In a case, if it is a bounded sequence (in particular, a finite 
sequence), then we add the last term to be +oo. Consider rectangles 

Tij = {{s, t) £ M 2 : U < s < t i+1 , tj<t< t j+1 } . 

Denote G{6) = {(t,y) : y = 9{t)}. 

By the construction, the rectangles have the following properties: 

- The set of all rectangles is at most a countable set; 

- The intersection G{6) n Tj is empty or contains a monotone part of the 
graph G{6). 

If G{9) H Tij is empty then we say Tj is empty. 

Now we shall construct the set 7baric(#)- Fix i, j such that the rectangle 
Tj is not-empty (an empty rectangle does not give any contribution to the 
set 7baric(0))> then we have 

Tba,nc(9)^7kj = a curve giving an one-to-one corespondence between [tj,tj + i] 
and [t k ,t k+1 ] if T ik / 0, k = 1, . . . , j - 1. 

Thus we have 

7baric(#) = U (Tbaric(^) H 7^) . 
kj 

Since there are a countable set of rectangles and in each rectangle we may 
have at most a curve which has Lebesgue measure zero (because, these curves 
give one-to-one correspondences), we conclude that the set 7b a ric(#) a ls° has 
zero Lebesgue measure. □ 

Example 4'. Consider the CEA E^ s '^ constructed in Example 4 by a 
family of invertible lower (or upper) triangular matrices t > 0. 

Theorem 4.9. For any pair of time (s, t) the n-dimensional evolution al- 
gebra E^^, constructed by a family of (lower or upper) triangular invertible 
matrices is baric. Moreover, E^ s ^ has a weight function a(x) = Ain$x n , 
where i = 1, ... ,n are diagonal entries of M\ s,t ^ = A^^lM) -1 . 
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Proof. It is known that the standard operations on triangular matrices con- 
veniently preserve the triangular form: the sum and product of two lower 
triangular matrices is again lower triangular. The inverse of a lower tri- 
angular matrix is also lower triangular, and of course we can multiply a 
lower triangular matrix by a constant and it will still be lower triangular. 
This means that the lower triangular matrices form a subalgebra of the ring 
of square matrices for any given size. The analogous result holds for upper 
triangular matrices. Using these properties we get that M\ s,t ^ is also a trian- 
gular matrix. Moreover, since A™ is invertible, its determinant is non-zero 
for all t. Thus 

n 

det(M [s > t] ) = \{M^ t] = det(^ H )det((^ M )- 1 ) ^ 0. 

4=1 

Consequently, all diagonal entries of the matrix are non-zero. In particular, 
M^nn 7^ 0, and Theorem 13.21 completes the proof. □ 

Corollary 4.10. The CEA E^ s,t ^ constructed by triangular invertible ma- 
trices has not baric property transition. 

4.11. Absolute nilpotent elements transition. The element x of an 
algebra A is called an absolute nilpotent if x 2 = 0. 

Let E = W n be an evolution algebra over the field R with structural 
constant coefficients matrix Ai = (a^), then for arbitrary x = ^2 i Xiei and 

V = Ui e i G ^™ we nave 



xy = X] ( X] a *i Xi ^ I e i' x2 = ^ ( X] a 

j \ i / j \ i 




J 

For a n-dimensional evolution algebra M. n consider operator V : 
x h-> V(x) = x' defined as 

n 

(4.2) x'j = ^2 a v x h 3 = !> • • • > n - 

i=i 

This operator is called evolution operator [11 1. 

We have V(x) = x 2 , hence the equation V(x) = x 2 = is given by the 
following system 

(4.3) ^2 Oijx 2 = 0, j = 1, . . . , n. 

i 

If det(A^) 7^ then the system f|4.3j) has unique solution (0, ...,0). If 
det(At) = and rank(A^) = r then we can assume that the first r rows of 
M. are linearly independent, consequently, the system of equations (|4.3|) can 
be written as 

n 

(4.4) x 2 = - ^2 d ij x % i = l,...,r, 

j=r+l 
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where = 'jffff with M r = K') ij=lv .., r , 



( a\\ ... aj-1,1 dji ... a r i ^ 

ai2 • • • Oj-1,2 %'2 Oj+1,2 • • • «r2 



\ ^lr • • • 0<i— \ % t Q'jr ^i+l,r • • • Ct rr J 

An interesting problem is to find a necessary and sufficient condition on 
matrix D = (dij) i=i,...,r under which the system (|4.4p has unique solution. 

j'=r+l,...,n 

The difficulty of the problem depends on rank r, here we shall consider the 
case r = n — 1. 

Proposition 4.12. 1) If det(A4) 7^ then the finite dimensional evolution 
algebra W 1 has unique absolute nilpotent (0, ...,0). 

2) If det(.M) = and rank(M) =n—l then the evolution algebra R n has 
unique absolute nilpotent (0, ...,0) if and only if 

(4.5) det(M ion ) • det(M n -i) > 0, 

for some io € {1, . . . ,n — 1}. 

Proof. 1) Straightforward. 

2) If rank(.M) = n — 1 then from (|4,4p we get 

From (14. 6p it follows that the condition (I4.5P is necessary and sufficient to 
have unique solution (0, . . . , 0). □ 

For a CEA E^ with matrix M\ s ^ denote 

Tnii = {(s,t)£T:E^ has unique absolute nilpotent}, 7^ = 7~\ 7^a- 

The following theorem gives an answer on problem of existence of "unique- 
ness of absolute nilpotent element" property transition. 

Theorem 4.13. 1) There are CEAs which have not "uniqueness of absolute 
nilpotent element" property transition. 

2) There is CEA which has "uniqueness of absolute nilpotent element" 
property transition. 

Proof. Denote d(s,t) = de^Al^). By equation (pT2|) we get 

(4.7) d(s,t) = d(s,T)d(r,t), for all r, s<r<t. 

As it was mentioned above, the equation (|4.7|) is known as Cantor's second 
equation. 

1) The equation (|4.7p has solutions d(s,t) = ^fj, where &(t) / is an 
arbitrary function. Thus for such solutions we conclude that if d(so,to) ^ 
for some (sq,£o) then d(s,t) ^ for any (s,t). Consequently, corresponding 
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CEAs have not "uniqueness of absolute nilpotent element" property transi- 
tion. 

2) Note that the equation (|4.7j) has solution d(s, t) = f(t), where f(t) = l 
for t < 1 and f(t) = otherwise. For this solution we have d(s,t) = l,s < 
t < 1 and d(s,t) = 0, t > 1. For some t > 1 one can construct a matrix 
.M^'*! which does not satisfy uniqueness condition mentioned in part 2) of 
Proposition 14.121 Indeed let us consider the matrix M\ s ^ = ^o|*-'^ 

with entries as in (12. 4ft . The second equation of the system (12.5P has a 
solution: 



i,i=l,2 



1, if s < t < 1 
0, if t > 1. 



/3(a,i) 

Substituting this solution in the first equation of (I2.5P we obtain 

a(s,i) = 



i(j(t) -ip(s), if s < t < 1 
$(*), if t > 1, 



where ?/> and g are arbitrary functions. The corresponding matrix has the 
following form 

= =-[ ), if s < t < 1, 



2 



m-^(s) 2-^(t)+v(s) 



and 



(4.8) M^ = -( , if t>l. 



We have 



d(s,t) = det(M 



[*,tb 



1, if s < t < 1 
0, if t > 1. 



Assume #(£) 7^—1 then for (|4. 8j) the equation (|4.6p has the form 

2 l-g(*) _2 

1 ! + <?(*) 

This equation has infinitely many solutions if |<?(i)| > 1 for some t > 1. 
Thus corresponding CEA has "uniqueness of absolute nilpotent element" 
property transition. □ 
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Now let us construct the set T n u for Examples 1-5: It is easy to see that 



det(M [s ' t] ) 



e~ 3At , for Example 1 

(A//)*, for Example 2 

for Example 3 

niLi-Mif 1 , for Example 4 
1, for Example 5 



Thus in each one of the considered examples we have det(Ai^ s '^) ^ 0, con- 
sequently, Tnu = T, i.e. the CEAs constructed in Examples 1-5 have not 
"uniqueness of the absolute nilpotent element" property transition. 

There are CEAs which have infinitely many absolute nilpotent elements 
independently on time. For example, take M\ s ^ with identical rows 
(^|y, 0, 0, . . . , 0), where <3? is an arbitrary function with <£(i) ^ for all t. 
It is easy to see that this matrix satisfies the equation (j2.2j) . hence it deter- 
mines a CEA, E^ s,t \ which has infinitely many absolute nilpotent elements: 
(0, X2, ■ ■ ■ , x n ), where X2, ■ ■ ■ ,x n G M are arbitrary numbers. Thus for this 
example we have 7~ n u = 0, T® u = T '■ In other words the CEA has not 
"non-uniqueness of absolute nilpotent element" property transition. 

Remark 4.14. These examples (Examples 1-4) of "uniqueness of nilpotent 
element" property transition of CEAs with time-parameter are similar to 
the "uniqueness of Gibbs phase" property transition, i.e. phase transition 
of physical systems with respect to temperature-parameter, T > 0. Usually 
there is a phase transition if the temperature is very low (T ~ 0) or if it 
is very high (T ~ +co) (see f^). Example 5 is an analogue of a physical 
system which has unique (Gibbs) phase for any temperature. There a lot of 
examples of such physical systems (see e.g. \B$)- 

4.15. Idempotent elements transition. A element x of an algebra A is 
called idempotent if x 2 = x; such points of an evolution algebra are especially 
important, because they are the fixed points (i.e. V(x) = x) of the evolution 
operator V, ()4.2j) . We denote by Zd(E) the idempotent elements of an 
algebra E. Using (|4.2p the equation x 2 = x can be written as 



(4.9) x j = ^ j a ij x 2 , j = l 



The general analysis of the solutions of the system (j4.9j) is very difficult. 
We shall solve this problem for the CEA t > 0, corresponding to the 
Example 2. In case of Example 2 the system (14. 9p has the following form 



(4.10) 

where A > 0, /j, > and t > 



2x = (A* + n l )x 2 + (A* - ^)y 2 ; 
2y = (A* - fi f )x 2 + (A* + ^ l )y 2 , 
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Case A = \i. It is easy to see that if A = \i then the system (14. 1QH has only 
four solutions = (0,0), z\ = z x {t) = (0,A~*),2 2 = z 2 (t) = (A _ *,0),z 3 = 
z 3 (t) = (\- t ,\- t ). 

Case A^ju. For A//i the solutions and z 3 still exist. If x = or y = 
there is no any new solution. Thus we consider the case xy ^ 0. Denote 

x A* -ax' (A//x) f - 1 

U y' 7l J A* + ^* (A/^ + l' 

For £ > it is easy to see that if A < fx then —1 < 7(t) < and if A > /x 
then < j(t) < 1. Note that for t = there is no any new solution. From 
system (|4,10p we get 

(4.11) 7 (i> 3 - u 2 + u - 7 (t) = (u - 1) {j(t)u 2 + ( 7 (t) - l)u + 7 (i)) = 0. 

Subcase A < /x. In this case for any i > the equation (|4.1ip has three 
solutions 



(A 19 , , l-7(0±v / l-2 7 (t)- 37 2 (t) 
(4.12) „ 1 = 1, u ± = ^ . 

Subcase A > fi. In this case the number of solutions to the equation (j4.11[) 
varies by 7, i.e. 



(4.13) solutions to (liTTT) 



if i < 7 (t) < 1; 
,«+ if < 7 (t) < ±, 

where ti± are defined in (I4.12|) . 

Now we shall describe x, y corresponding to the solutions of (I4.1ip . The 
case u = 1, i.e. x = y does not give any new solution. For u = u± we 
have x = u±y, substituting this in the second equation of (|4.10p after simple 
calculations we get the following two non-zero solutions to (|4.10p : 



// ± JX t {2n t - A*) A* - /i* 

s ± = -T77T7 TtttS 5 V± ~ 



At* (A* ± V / A*(2ix* - A')) p*(A* ± -v/A*(2/x* - A*)) 

Note that cc±, y-j- are well defined for any A / /j. For A > \i we have critical 
time 

(4.14) t c 11,2 



In A — In n 

which is the unique solution to the equation 7 (i) = ^. 
Thus we have proved the following 

Proposition 4.16. We have 

'{0,21,22, 23}, */ A = Ml 

[*]) = < {°' z 3, (a;-,y-), (^+,2/+)}, «/ A</i; 

{0,z 3 }, i/ A > a*; * > t c 

k {0,2 3 , (x_,y_), if A > /i, i < t c . 
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This proposition gives a very interesting "a fixed set of idempotent ele- 
ments" property transition, i.e. we have 

Corollary 4.17. The CEA constructed in Example 2 has not "a fixed 
set of idempotent elements" property transition if A < \i; it has such property 
transition if X > [i. Moreover the transition point (the critical time) is t = t c 
defined by formula l[4.14\ )- 

Remark 4.18. There are exactly solvable models in statistical mechanics, 
here an imprecise notion of "exactly solvable" as meaning: "The solutions 
can be expressed explicitly in terms of some previously known functions" is 
also sometimes used JI]/. In such models, for example, the critical temper- 
ature can be expressed explicitly. Comparing this with our examples of a 
property transition we also can say that a property transition of a CEA is 
exactly solvable if the critical time can be found exactly. Thus our Example 
2 is exactly solvable for investigation of properties of idempotent elements. 
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